Model-independent test of gravity with a network of ground-based gravitational-wave 

detectors 
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The observation of gravitational waves with a global network of interferometric detectors such 
as advanced LIGO, advanced Virgo, and KAGRA will make it possible to probe into the nature 
of space-time structure. Besides Einstein's general theory of relativity, there are several theories of 
gravitation that passed experimental tests so far. The gravitational-wave observation provides a new 
experimental test of alternative theories of gravity because a gravitational wave may have at most 
six independent modes of polarization, of which properties and number of modes are dependent 
on theories of gravity. This paper proposes a method to reconstruct the independent modes of 
polarization in time-series data of an advanced detector network. Since the method does not rely 
on any specific model, it gives model-independent test of alternative theories of gravity. 

PACS numbers: 04.80.Cc,04.30.-w,04.80.Nn,07.05.Kf 



I. INTRODUCTION 

In recent years, direct detection experiments of a grav- 
itational wave (GW) have been well developed and the 
first generation of a kilometer-scale ground-based laser- 
interferometric GW detector has accomplished its design 
sensitivity. Although the first detection of the GW has 
not been achieved yet, the null detection has yielded sci- 
entific results [iJ-IH . The next-generation interferometers 
such as advanced LIGO @, advanced VIRGO Q, and 
KAGRA Q will be in operation in the coming five years 
and will bring valuable information about astronomical 
compact objects. 

The direct observation of the GWs will also provide 
a unique opportunity to test the theory of general rela- 
tivity (GR), through the propagation speed, waveforms, 
and polarization modes of GWs. In GR, a GW has two 
polarization modes (plus and cross modes), while in a 
general metric theory of gravitation, the GW is allowed 
to have at most six polarizations 0, [l(|. In modified 
gravity theories such as the scalar-tensor theory [ill, E3 
and f(R) gravity [U, EH, additional scalar polarizations 
appear (For more rigorous treatment of the polarizations 
with the Newman-Pcnrosc formalism, see Imllfi). On 
the other hand, in bimetric gravity theory |17j and mas- 
sive gravity theory [IH, [lj| , there appear at most six and 
five polarization modes, respectively, including scalar and 
vector modes 0, [2(| • If the additional polarizations are 
found, it indicates that the theory of gravitation should 
be extended beyond GR and excludes some theoretical 
models, depending on which polarization modes are de- 
tected. Thus, the observation of the GW polarizations is 
a powerful tool to probe the extended law of gravity. 
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Currently, there are few observational constraints on 
the additional polarization modes of GWs. For the 
scalar GWs, the observed orbital-period derivative of 
PSR B1913+16 agrees well with predicted values of GR, 
conservatively, at a level of 1 % error [2l| , indicating that 
the contribution of scalar GWs to the energy loss is less 
than 1 %. However, it is important to cross-check the 
existence of the number of propagating degree of free- 
dom directly by GW detection experiments, since these 
experiments probe in a weaker regime of gravity than 
a binary pulsar and at different distance scale. So far 
some authors have been proposed the method for sepa- 
rating a mixture of the polarization modes of the GW 
background and detecting non-Einsteinian polarization 
modes with pulsar timing array [22j and with ground- 
based laser- interferometric GW detectors [23|, [24|. To 
decompose and reconstruct the polarization modes, the 
number of the independ! cnt signals of detectors should 
be more than that of polarization modes. In the above 
method with GW detectors, the number of polarizations 
are assumed to be three, since they considered a stochas- 
tic background. 

Now we expect that more GW detectors will be in 
operation in the future, in total five detectors including 
two advanced LIGOs, advanced VIRGO, KAGRA, and 
IndlGO [25j]. Then in principle we can reconstruct five 
polarization modes (all modes that GW detectors can 
separately detect) of a chirp and burst GW. In this pa- 
per, we propose a method to separate and reconstruct 
the polarization modes using observation data of multi- 
ple interferometric GW detectors. The proposed method 
does not need any theoretical waveform of the polar- 
ization modes of a GW, and therefore can be a model- 
independent probe of testing various alternative theories 
of gravity. 

The organization of this paper is as follows. In sec- 
tion [IT] we first describe antenna pattern functions of the 
tensorial and the scalar polarization modes, and then 
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construct the algorithm to separate and reconstruct the 
polarization modes based on the coherent network anal- 
ysis. In section IIV A| we describe how the algorithm is 
implemented and show an example of reconstruction of 
the polarization modes. In section IIVI we made simu- 
lations of the reconstruction of the scalar polarization 
mode using simulated data of LIGO Hanford, LIGO Liv- 
ingston, VIRGO, KAGRA. We devote the last section [V] 
to the summary of this paper. 



II. ANALYTICAL METHOD OF 
RECONSTRUCTION 

In this section, we provide the method separately de- 
tecting and reconstructing more than three polarization 
modes of a GW, which often appear alternative theories 
of gravity, with a coherent network of ground-based de- 
tectors. 



A. Polarization modes of a gravitational wave 

In general, a metric gravity theory in four dimensions 
allows at most six polarization modes of a GW 0, ■ 
Let us define a wave orthonormal coordinate that are 
constructed by a unit vector k directed to the propaga- 
tion direction of a GW and two unit vectors eg and 
orthogonal to k and each other. With these vectors, the 
polarization modes are defined as 



e + = eg ® eg - ® , 


(1) 


e x = eg <g) e + e <E> eg , 


(2) 


e° = eg ® eg + e^ <g) , 


(3) 


e e = V2 k <g> k , 


(4) 


e x = eg ® k + k ® eg , 


(5) 


e lJ = e^ ® k + k ® , 


(6) 



where the symbol <£> denotes a tensor product. The +, x , 
o, t, x, and y polarization modes are called plus, cross, 
breathing, longitudinal, vector-x, and vector-y modes, re- 
spectively. According to rotation symmetry around the 
propagation axis of the GW, the + and x modes are iden- 
tified with tensor-type (spin-2) GWs, the x and y modes 
are vector-type (spin-1) GWs, and the o and I modes 
are scalar-type (spin-0) GWs. For the more detailed in- 
troduction, see [23. With the sky direction of a GW 
source O = —k, then a GW with the six polarizations is 
expressed as 

h ij (t,Cl) = J2h A (t)et j (n), (7) 

A 

where A — +, x, o, £, x, y. 



B. Antenna pattern functions of polarization 
modes 

Antenna pattern functions of the scalar and vector 
modes have been derived in [2^, [2(| [27j . We briefly sum- 
marize the results here. 

Suppose m intcrferomctric GW detectors are in oper- 
ation. The GW signal of the 7-th detector is written as 

Z I (t,£l)=J2Fi A (fyh A (t) , (8) 

A 

where Fj(Cl) is the antenna pattern function of the 7-th 
detector defined as 

F I A (fl) = ef j (Q)4 j . (9) 

Here, di is a detector tensor defined as 

di := - [uj <g> ui — vi vi] , (10) 

where uj, vi are unit vectors along with arms of the I-th 
detector. In a spherical coordinate (6i,4>i) fixed to the 
7-th detector, since the detector coordinate (u,v,w) and 
the wave coordinate (eg, e,p, 0) are related by 

eg = iii cosOi cos 4>i + vi cos9i sin</>/ — wi sin0/, (11) 
= — ui sin (f>i + ii cos (pi , (12) 
Cl = ui sin9i cos 4>i + iii sin#/ sin </>/ + ibi cos 61 , (13) 

the angular pattern functions for each polarization are 

F+(0) = -(1 + cos 2 6»/) cos 20/ , (14) 
Ff(£l) = -cos6»/sin20/ , (15) 
Fi(Cl) = -i sin 2 6»/ cos 20/ , (16) 

Ff(h) = ^ sin 2 6>/ cos 20/ , (17) 

F?(£l) = -i sin 26*/ cos 20/ , (18) 

F?(Cl) = sin 6»/ sin 20/ . (19) 

Note that the wave coordinate has a rotational degree of 
freedom tp about the axis. The general formulas keep- 
ing "0 arbitrary are provided in [23j. From the expres- 
sions, one can see explicit rotational symmetries about 
ip for the scalar, vector, and tensor polarization modes. 
In this paper we take "0 = without loss of generality. 
Also we note that the above expressions are valid when 
the arm length of the detector, L, is much smaller than 
the wavelength of observed GWs, A g , i.e., L<A g . This 
condition is well satisfied for ground-based detectors we 
consider in this paper. The angul ar p attern functions for 
L > Aff have been derived in |28|, [29| (For pulsar timing, 
see [30, [3lJ])- Figure [1] shows the scalar, vector, and ten- 
sor antenna patterns as a function of the longitude and 
the latitude. 
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FIG. 1: The left top plot is an antenna pattern of LIGO interferometric detector at Hanford to a scalar GW (F 2 ) 1 ^ 2 . The 
X-axis is the longitude and the F-axis is the latitude. The right top plot is an antenna pattern function of the detector to the 
tensorial modes (F 2 + F 2 ) 1 / 2 . The bottom plot is an antenna pattern function of the detector to the vector modes (F 2 + F 2 ) 1 ^ 2 



The angular pattern functions for scalar modes in 
Eqs. (|16p and (fTT)) are degenerated. Thus, these GW 
signals of the scalar modes cannot be distinguished in 
GW observation and what we actually detect is the sig- 
nal combination ^Fj{h (t) — \/2he(t)}. In addition, the 
scalar logitudinal mode often appears as a massive ghost 
mode in a modified gravity theory and should be removed 
from the theory From these reasons, hereafter we as- 
sume that the longitudinal mode is absent. In addition, 
we assume that all polarization modes are massless since 
separate analysis is needed for the mode that propagates 
with the speed much different from that of light. There- 
fore, we consider in total five massless polarization modes 
(+, x, o, x, y) of a GW in the following sections. 



C. Coherent network analysis 

The coherent network analysis is an algorithm to find 
a GW signal in the data b y c ombining all available de- 
tector data coherently (|32|-|37| and therein). The detec- 
tor output is a combination of GWs in the polarization 
modes weighed by the pattern function of each polariza- 
tion mode. In the coherent network analysis, the sky 
location of the GW and the waveforms in all polariza- 
tion modes are reconstructed by inverting the set of the 
detector responses. 

Data xi(t) from the 7-th detector is 



(20) 



where £r(i) and r)i(t) are the GW signal and noise of 
the 7-th detector. The noise is assumed to be Gaussian 
distributed. The arrival time of a GW at each detector 
is delayed depending on the geographical locations of the 
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GW detectors. If the relative time delay with respect to 
a reference time to taken at the center of the Earth is 
defined as tj(c/),9), the arrival time can be redefined as 
t = t + T I (<l>,6). 

All the detectors being taken into account, Eq. (f20|) is 
written as 



where 



x = Fh + r\. 



X — ' ' ' ) Xrn) 7 

/ F+ F* F° Ff Ff \ 
F = : : : : : 

jp+ px po pi jpy 

\ m m m m m / 

:= (F + ,F x ,F ,F x ,F y ), 
h = (h + ,h x ,h ,h x ,h y ) T , 



and 



XI 

Fj A 
h A 

m 



(xi(fo), ■ ■ ■ ,£/(/jv-i)) 



N-l 



(M/o), 

(m(fo), 



,h A (fN-i)) 

,fjl(fN-l)) ■ 



(21) 



(22) 



(23) 
(24) 
(25) 



(26) 
(27) 
(28) 
(29) 



The subscripts run as / = 1 , • • • , m and A = 
+, x, o, x, y. The components are defined in Fourier 
space as xjjfj) := 5i(/ 3 -)/VS2(/j)> Ff&fj ) ■= 
F^VVW), and mifj) := mUiVy/W^, J = 
0, • • • , N — 1. i(/j) is defined as j-th component of the 
Fourier transform of x(t): 



N-l 



(30) 



W is the number of samplings, A* is a sampling period, 
Af := (NAf)^ 1 is a frequency resolution, and fj = jAf. 
xi(fj), the whitened x(fj), is obtained by dividing by 
the power spectrum density S^(fj) of 7-th detector at 
a frequency fj. There are several purposes to whiten 
detector data. The sensitivity of an interferometric GW 
detector is limited in a freq uency- dependent way through 
a diverse noise budget [38|. Since multiple data streams 
have different transfer functions in the Fourier domain, 
one needs to whiten the data streams in order to combine 
them coherently Otherwise insensitive frequency regions 
or insensitive detectors in a detector network dominate 
when the detection statistics is calculated. If data is 
whitened, these frequency regions or detectors are less- 
weighted naturally. Another important purpose of the 
data conditioning is to make the data un-correlated be- 
tween samples by a whitening filter, and to mitigate the 
effect of instrumental artifacts in detector data. For in- 
stance, noise artifacts appear in multiple detectors with 



the same frequency regions, can be correlated noise be- 
tween detectors and this makes the detection efficiency 
reduced. 

We first consider the reconstruction of waveforms in an 
ideal case without noise. In general, F is not a squared 
matrix with full rank, hence we introduce the Moorc- 
Pcnrose pseudo-inverse matrix M as 



Mh = F T x, where M := F T F. 



(31) 



If the detectors are not all co-aligned, M is an invertible 
A x A matrix. Multiplying the equation by the inverse 
of M . we get 

h = F^x, where F f := M^ 1 F T . (32) 

The inverse matrix M~ l can be found from the formula 

f 



M = 



dct(M) 



adj(M), 



(33) 



where adj(Af) is the adjoint matrix and det(Af) is the 
determinant of M. If (i,j)-th cofactor is defined as 
Caa' = \^}{M)\aa' i we finally obtain 



h A 



H a x , 
1 

dct(M) 



AA' 



(34) 
(35) 



As discussed in [2J], the factor det(iW) of this formula 
plays an important role in separating the polarization 
modes. If there is degeneracy in the antenna pattern 
functions of a detector network and det(M) = 0, we 
cannot reconstruct the polarization modes at all. For- 
tunately, this is not the case for current ground-based 
detectors. 

We use a maximum likelihood method to estimate h 
from the data. The maximum likelihood method maxi- 
mizes 



L[h) : = -\\x-Fh || 2 , 

= - || x- FF^x || 2 , 

where || ■ || is defined by 



m N-l 

££M^/)I 2 ^/ 

1=1 3=0 



(36) 
(37) 



1/2 



(38) 



Introducing Q 
obtain 



I - FF^ and using Q in Eq. (|37p . we 



L[h] 



Qx 



(39) 



We note QF = 0, which means that Q projects onto the 
null space of FF T . From Eq. (|39p . one can see L[h] is 
equivalent to the null stream energy [39j . Suppose the 
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true source location is 17, 



L[h] = - || Q(F(n s )h + Tj) 



- || QF(Cl s )h || 2 -2Re 



j 

(40) 



Here, we used the relation Q 2 = Q and denoted fre- 
quency components as {rfQrj)j = [rf QT]](fj) . As 
Ref. [35j shows, if we take the expectation value of 
the likelihood function in a stationary noise case, the 
first term in Eq. (|40|) is the same as the current value 
|| QF(Ct s )h || 2 and the second term vanishes. The third 
term at j-th frequency bin is 

m 

i,j 



i.j 



(41) 



where Su is Kronecker delta. Since the noise in the detec- 
tors is whitened, the noise can be regarded as the Gaus- 
sian white noise with its variance a 2 := Y]j a 2 Af, where 



07 is 



defined by {{ffjYffj) 



SijcTj- Wc note a 2 has the 
same value for all j-th frequency bin since the noise ij 
is whitened. Transforming Q J to the diagonal form, we 
find tr(Q J ) = m — n and finally obtain 



(L[h}) 



QF(Cl s )h || 2 -{m-n)<j 2 



(42) 



where n is the total number of polarization modes. Since 
QF(Ct s ) = at the true location, L[h] become maxi- 
mum. 

We emphasize that the algorithm does not specify any 
alternative theories of gravity and can be applied to any 
number of polarization modes if it is less than the num- 
ber of detectors. Therefore, this approach is a model- 
independent probe alternative theories of gravity. 



III. APPLICATION TO THREE 
POLARIZATION MODES 

A. Scalar GW in Brans-Dicke theory 

There are several simulations on the spherically sym- 
metric core collapse in the Brans-Dicke theory , which 
is a class of scalar-tensor theory. The scalar gravita- 
tional waveform is consistent with each simulation 0- 
[45[ . Therefore we use a waveform simulated by Shibata 
et al. [40:| . In their simulation, they assumed the Brans- 
Dicke theory with the scalar field coupling to gravity, 
^bd = 50 and 500, and calculated scalar gravitational 
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FIG. 2: A scalar gravitational waveform simulated by 
Shibata et al. [4p(| with the mass of the progenitor IOMq, 
wbd = 40000, the distance from the earth lOMpc. The unit 
of the time t is 4.93 x 10~ 5 (A//10A/ o ). $ is the normalized 
strain amplitude. 



waveforms. They found the scalar fields are linearly 
scaled with cjbd , and the scalar gravitational waveforms 
with different cjbd are almost the same by scaling wbd- 
We therefore extrapolate the waveforms to larger values 
of wbd • The simulated scalar GW signal h is 



h(t) = 1.25 x 10" 



M 



flOkpcX 



ioM Q ; v R J V wbd ) 



(43) 

where M is the mass of a progenitor, R is the distance 
from the earth. The shape of the waveform $(<) is in 
Figure[2j Since the unit of the time t in Figure[2]is 4.93 x 
1O~ 5 (M/1OM0), The duration of the signal is linearly 
dependent of a mass of a progenitor. 



B. Sensitivity to the scalar mode 



Sensitivity of the designed Japanese interferometric 
GW detector KAGRA jifj to a scalar GW simulated by 
Shibata et al. [4(J is estimated. We should note that more 
realistic estimation of the detectability needs more real- 
istic simulation. However, we use Shibata's result for the 
purpose of demonstration. The top plot in Fig. [3] shows 
h rss -f c of GWs from spherically symmetric core collapse 
supernovae located at lOkpc from the earth with wbd — 
500, 1000, 2000, 4000, 8000, 16000, 40000, 80000, 160000. 
h TSS is the root sum square of h defined as 




1/2 



(44) 
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(SNR) as a function of cjbd- The SNR is defined as 



SNR 



Sn(f) 1 



1/2 



(47) 



We note that the SNR defined here is that of the matched 
filter method [48| , which is the optimal case with respect 
to the antenna pattern function of the detector. The 
SNR is estimated by 



SNR 



VfcSn(fc) ' 



(48) 



where h c is the characteristic strain amplitude [47 



h r := 3 



Sn(fc) 
' Sn(f) 



1/2 



(49) 



FIG. 3: Blue open circles are h vss -fc of spheri- 
cally symmetric core collapse supernovae at the distance 
lOkpc with the progenitor mass 20Mq and cjbd = 
500, 1000, 2000, 4000, 8000, 16000, 40000, 80000, 160000. h ISS , 
f c are defined in the text. Magenta open circles are for the 
progenitor mass 5Mq . The bottom plot is ojbd-SNR. SNR is 
defined in the text. The red line is the one with the mass of 
the progenitor 20Mq, the magenta line 5M©. 



A = 



x, o, and f c is the characteristic frequency [4 

fc- = 



E.M/)W) /d/ x (45) 



Sn(f) 

EA~hA(f)h* A (f) 

Sn(f) 



df 



(46) 



where x(f), S n (f), h A (f) are detector output, the noise 
power spectrum density, the plus, cross, and scalar modes 
of a GW in Fourier domain. h* A expresses its complex 
conjugate. Note that here we are considering GWs from 
spherically symmetric core collapse supernovae, which 
does not radiate tensor modes, namely, h+(t) = hx(t) = 
0. h rss of the supernovae with the progenitor mass of 
20Mq at the distance lOkpc are shown with blue open 
circles. The magenta open circles are the same except 
the mass of the progenitor is 5M . These plots are on 
the current design sensitivity curve of KAGRA. The bot- 
tom plot in Fig. |3] shows estimated signal-to-noise ratio 



The red and magenta lines are the ones with the mass of 
the progenitor 20M Q and 5M Q , respectively. The current 
upper limit on wsnis > 40000 by Cassini's observation in 
the solar system [49|]. If a scalar GW with ojbd = 40000 
is radiated from a spherically symmetric core collapse 
supernova at 10 kpc from the earth, KAGRA can detect 
the signal with SNR of 35 for the mass of 20M Q progen- 
itor, 7 for the mass of 5M progenitor. In order to set 
upper limit on ubd, one has to know the distance of a 
source because the distance and wbd influence the ampli- 
tude of the GW at a detector in the same way. Therefore 
the two parameters are degenerated. It is important to 
combine usual astronomical observation such as electro- 
magnetic observations and neutrino observations with 
the GW observation, so called the multi-messenger ob- 
servation. Suppose it find a supernova occurs within our 
Galaxy through the multi-messenger observation. The 
above results indicates one may set stronger upper limit 
on wbd, or if the scalar GW is detected, one can deter- 
mine wbd- 



C. Reconstruction of three polarization modes 



In the case the tensor and scalar modes exist, h + , 
h x , h are reconstructed as the linear combination of 
the antenna patterns of all the three polarizations. H A 
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(A = +, x, o) can be explicitly calculated from Eq. ([35 
1 
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+ 



-[-(F+ xF )-(F x xF )F + 



det(M)^x X ^-( F x X ^ F + 
(F x x F ) ■ (F + x F )F X 
(Fx xF ).(F + xF x )F ], 
1 

det(M) 
(F + x F ) • (F + x F )F X 
(F + xF )-(F + xF x )F ], 

dWM) [(F+XFx) ' (FxXF ° )F+ 
(F + xF x ).(F + xF )F x 

(f + xF x ).(F + xf x )f ]. 



(50) 



(51) 



(52) 



The symbols x and • denote the outer product and the in- 
ner product, respectively. We note if det(Af) = 0, which 
happens when detectors are co-aligned, the reconstruc- 
tion of any polarization mode fails. However, this is not 
the case for the current ground-based detector network. 



IV. SIMULATION OF RECONSTRUCTION OF 
POLARIZATION MODES 

A. Implementation 

For the data conditioning, the linear predictor error 
filter is used, assuming data is a stationary stochastic 
process and can be expressed by an autoregressive model 
with order of P. If predicted data stream can be written 

as 



x'(sA t 



p 

£ 



c(pA t )x((s-p)A t ) 



(53) 



the linear predictor error filter finds parameters c(pA t ) 
which minimize the prediction error 



x w (sA t ) := x(sA t )-x'(sA t ) 



p 

Y,b(pA t )x((s-p)A t ) 

p=0 



(54) 
(55) 



using Yule-Walker method [5C§. Here b{pA t ) := S p q — 
c(pA t ), c(0) := 0. This filter has a group delay of the 
phase especially near narrowband spectral feature and 
the timing error between the whitened data of different 
detectors is introduced. To cancel the timing error, the 
linear predictor filter is first applied causally and then 
anti-causally, 



x w (f) = B(f)B*(f)x(f), 



(56) 



where x w (f) is the discrete time Fourier transform of 
whitened data, 



B(f) 



OO 

£ 

m=0 



6(pA t )exp(-i27r/p//.), 



(57) 
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FIG. 4: The sensitivity curves. The i-axis is frequency and 
the y-axis is amplitude spectrum density. Detectors consid- 
ered are H, L, V, K with the same sensitivity which is similar 
to the design sensitivity of advanced LIGO. 



f a is a sampling frequency of the data [5l|. The filter 
coefficients b(pA t ) is estimated using a stationary data 
segment. For this purpose, we use a few seconds data seg- 
ment that docs not have non-stationary noise in the im- 
plementation. The data segment will not used for search 
for GWs. We than construct the finite impulse response 
(FIR) filter with the coefficients b{pA t ) and pass the data 
through the FIR filter to whiten it. 

Set of the whitened data x of the multiple detectors 
with equal length is passed to the coherent network anal- 
ysis pipeline. The value of the likelihood of the data 
which is maximized over all possible h + , h x , and h with 
durations less than or equal to the data length. The 
likelihood values are obtained function of 9 and d>. 
This two-dimensional output, S(0,<p), is called skymap. 
Finally reconstructed h +1 h x and h are evaluated the 
maximum point of the skymap. We note that the imple- 
mentation is based on the paper p36| . 



B. Reconstruction of the polarization modes 
without tensor mode signals 

We perform Monte Carlo simulations to reconstruct 
the scalar polarization of a GW. The network consisted 
of the 4 km LIGO Hanford (H), LIGO Livingston (L), 
VIRGO (V), and KAGRA (K) interferometers 1^ . For 
the detector noise amplitude spectral densities, we use 
the design sensitivity curves similar to the advanced 
LIGO detectors as given in [52| (see Fig. [3} and keep 
the locations and orientations the same as the real de- 
tectors. The stationary noise is generated 50 seconds by 
first using 4 independent realizations of Gaussian white 
noise and then passing them through FIR filters having 
transfer functions that approximately match the design 
curves. The generated data is sampled at 16384 Hz and 
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FIG. 5: Reconstructed polarization modes. The top plot 
is h+, the middle plot is h x , the bottom plot is h a - The 
x axis is time with the unit of millisecond since GPS time 
873651621. 973 14453 [s]. The black dashed line in each plot is 
the injected signal. 



then passed through the data conditioning pipeline. Be- 
sides downsampling the data to 2048 Hz by applying the 
same anti-aliasing filter to all data streams, the data con- 
ditioning pipeline applies time domain whitening filters 
that are trained on the first five seconds of data in which 
any injected signal is not included. In this simulation we 
put cut-off frequencies at 60 Hz and 400 Hz so that lower 
and higher frequencies are filtered. 

We assume that the GW source is a spherically sym- 
metric core collapse supernovae. The injected signals 
corresponds to a single source located at the right as- 
cension (RA) of 15 hours and the declination (DEC) 
of —60 degrees. The mass of the progenitor is 10M Q , 
^bd = 40000, the distance from the earth is lOkpc. we 
here assume that h x (t) and h+(t) are absent and that 
the simulated gravitational waveform is h . In the next 
subsection, we will perform the same simulation in the 
presence of h x (t) and h + (t). 

Fig. [5] shows the reconstruction of h+, h x , h , which is 
one of segments triggered by the pipeline. In the recon- 
structed time-series of h , the injected signal is clearly 
reconstructed. The injected h Q different from Figure [2] 
This difference comes from the fact that the low fre- 
quency region below 60 Hz and the high frequency region 
above 400 Hz in data are filtered in this simulation. 

In this paper, it is supposed that the gravitational 
waveform is known and the matched filter method is ap- 
plicable to the search. To evaluate the accuracy of the 
reconstruction, the SNR is calculated in the same way of 
Eq. gTJ): 



SNR 



SnA(f) 



1/2 




df 



(58) 



FIG. 6: The SNR of F h in H, L, V, and K, and recon- 
structed h as a function of wbd • 



where f)yt(/) is the reconstructed Ha in Fourier domain, 
fiA(f) is the injected Ka in Fourier domain, and S n A(f) 
is the power spectrum density of the reconstructed Ha- 
It should be noted that realistic searches for GWs from 
supernovae cannot use the waveform models because the 
models are still not well-established at present. The re- 
alistic search may use other indicators which do not use 
waveform models such as the excess power statistics [53j . 

The SNR of the reconstructed h which is band- filtered 
is 14.3. The SNR of h Q in H, L, V, K before the recon- 
struction is 4.9, 13.4, 16.4, 0.54. F is -0.26 for H, 0.42 
for L, —0.45 for V, —0.01 for K. The main contribution 
to the SNR is from V which has best antenna pattern 
among them. 

Fig.[B]shows the SNR of h in H (the magenta line), L 
(the blue line), V (the light blue line), and K (the green 
line), and reconstructed h (the red line) as a function 
of wbd- The reconstruction procedure were repeated 28 
times with different data set and the SNRs were calcu- 
lated at each trial and then they were averaged. Inter- 
estingly, the SNR of the reconstructed h is higher than 
the others. As shown in Eq. (|37[) . the maximum like- 
lihood method combines all detector data streams that 
are weighted by their sensitivities in the whitening pro- 
cedure. Multiple sensitive-detector data streams make 
contribution to reduce the variance of the reconstructed 
h and the SNR of the reconstructed h rises up. In re- 
gard to constraining wbd by practical GW observations, 
further studies including the difference of each detector 
sensitivity and the dependence of SNR on a sky position 
are needed. 



C. Reconstruction of the polarization modes with 
tensor mode signals 

In previous section we considered the reconstruction 
of the GW which contains only scalar mode. The al- 
gorithm can reconstruct all the modes if they exist. In 
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FIG. 7: Reconstructed h+, h x , and h a . The black dashed 
line on each plot is the injected waveform. 



this section we demonstrate the reconstruction of a GW 
which contains three modes, h +1 h x , h . Example of the 
reconstruction of the polarization modes h+, h Xl and 
h is shown in Fig. [7] The detectors H, L, V, K are 
assumed to have the same sensitivity as shown Fig. @] 
The set of detector noise is generated in the same way 
in section ITVl Regarding to injection signals, in this sim- 
ulation we will not use GW waveforms predicted by the 
Brans-Dicke theory because spherically symmetric core 
collapse supernovae do not produce tensor modes, and 
there is no realistic simulations on asymmetric core col- 
lapses in the Brans-Dicke theory. Instead, we will use the 
sine- Gaussian signals that are usually used for evaluation 
of GW burst searches. The injected signals correspond 
to a single source located at the right ascension (RA) of 
16.4 hours and the declination (DEC) of degrees. We 
assume that 



M*) 
h+(t) 

ho(t) 



AeM-(^f t) 2 /2Q 2 } sin(27r/ t), (59) 
Aexp[-(27r/ i) 2 /2Q 2 ] sin(27r/ i + tt/2), (60) 
Aexp[-(27r/ <) 2 /2Q 2 ] sin[27r/o(i + 0.003)], 

(61) 



where Q = 9 is the Q-value and /o = 235 Hz is the cen- 
tral frequency. The signal strength A is scaled so that 
the root-sum-square h ISS = 1.38 x 10~ 22 Hz~ 1/2 . The 
dashed lines in Fig. [7] are the injection signals. We eval- 
uate the signal strength by SNR. Table U shows the SNR 
of the reconstructed h + , h x , h Q , and the SNR before re- 
construction of the three polarization modes in H, L, V, 
K. 

If we define a mean SNR as 



SNR = 



SNR(h H 



SNR(h x ) 2 +SNR(/i 



(62) 





Reconstructed 


H 


L 


V 


K 


h+ 


23.3 


8.3 


9.6 


26.7 


24.5 


h x 


20.2 


15.4 


18.0 


15.4 


10.8 


h 


16.3 


13.3 


20.4 


10.0 


6.2 



TABLE I: The SNR of the reconstructed h+, h x , h a , and 
the SNR of the three polarization modes in H, L, V, K. The 
antenna pattern function (F + ,F X ,F ) is (0.25,-0.46,0.39) 
for H, (0.19,0.36,-0.41) for L, (-0.65,-0.38,-0.24) for V, 
(0.56,0.25,0.14) for K. 



even if several detectors in the detector network do not 
have sensitivities to certain polarization modes, sensitive 
detectors in the detector network can make improvement 
of the SNR. One indication to know the detector-network 
sensitivity to the polarization modes is the network an- 
tenna pattern skymap F T := v(X)jli F™+F™), Fs '■= 
TlE/LiO, F V := V(E?^F[ 2 + F^). Figured 
shows the network antenna pattern skymap of H-L-V-K. 
Comparing the antenna pattern skymaps of the polariza- 
tion modes for H (see Figure [JJ, for the antenna pattern 
skymaps of the tensor modes, the region where the value 
is above 0.5 is increased from 72.8% to 100%, for the 
vector modes, the region where the value is above 0.5 is 
increased from 69.1% to 100%, and for the scalar mode, 
the region where the value is above 0.25 is increased from 
31.4% to 94.2%. These shows the benefit of the use of 
multiple detectors in terms of the improvement of the 
SNR and the sky coverage. The latter means that angu- 
lar directions insentive to a certain polarization mode is 
removed and that the detector network can distinguish 
all polarizations. It should be noted that the matched 
filtering method on reconstructed polarization modes is 
proposed in (54| • This method benefit from the detector 
network directly. 



SNR is 20.1 for reconstructed modes, 12.7 for H, 16.7 for 
L, 18.7 for V, 15.9 for K. From the Table HI one can see 



V. SUMMARY 

This paper proposed a method to reconstruct at most 
five non-Einsteinian polarization modes, in addition to 
h + and h x , of GWs using a network of ground-based in- 
terferometric GW detectors. This method is applicable 
for testing alternative theories of gravity by searching for 
scalar and vector GWs. Since the method does not rely 
on any specific model of the proposed alternative theo- 
ries and does not need theoretical models of gravitational 
waveforms, the method can be a model-independent test 
of alternative theory of gravity. 

We overview the algorithm of the method. The de- 
tector responses is first whitened by the linear predictor 
error filter so that the processed data is sample to sample 
uncorrelated. Since the detector response to a GW is the 
linear combination of the polarization modes weighted 
with the antenna pattern functions, the reconstruction 
of the polarization modes is naturally formulated as an 
inverse problem. We solved the inverse problem by min- 
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FIG. 8: The left top plot is the network antenna pattern of H-L-V-K to a scalar gravitational wave. The X-axis is the longitude 
and the F-axis is the latitude. The right top plot is the network antenna pattern function to the tensorial modes. The bottom 
plot is the summed antenna pattern function to the vector modes 



imizing residuals by subtracting a reconstructed signal 
from the detector output data. 

We performed simulations of the reconstruction of the 
polarization modes as demonstrations. We used the sim- 
ulated gravitational waveform from a spherically sym- 
metric core collapse supernova in Brans-Dicke theory, 
which predicts only the scalar GW, and showed the scalar 
GW was well reconstructed. We also demonstrated that 
all gravitational waveforms were well reconstructed in 
the presence of two tensor modes in addition to a scalar 
mode, using sine-Gaussian waveforms. The mean SNR 
defined in Eq. (|62j) was calculated for the reconstructed 
polarization modes and the polarization modes in H, L, 
V, K before reconstruction. The mean SNR of the recon- 
structed polarizations were higher than the others, which 
means that sensitive detectors play an important role to 
reconstruct the waveforms and gain the SNR. For fur- 
ther study of the detectability, more realistic simulations, 
considering the sky position dependence of a source and 



other waveforms in a specific model of gravity theory, are 
encouraged. 

Finally we discuss the case that other polarization 
modes than the three polarization modes exist. In this 
case, because Q does not project onto the null space of 
FF T where F includes all polarization modes, L[h] has 
larger variance than expected noise variance in Eq. (|42[) 
even if (6, (f>) is the true source location. Therefore we 
can test the existence of other polarization modes. 
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